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ABSTRACT 

The distant type la supernovae data compiled by Perlmutter et al. (1999) are 
used to analyze the Cardassian expansion scenario, which was recently proposed 
by Freese and Lewis (2002) as an alternative to a cosmological constant (or more 
generally a dark energy component) in explaining the currently accelerating uni- 
verse. We show that the allowed intervals for n and z eq , the two parameters of 
the Cardassian model, will give rise to a universe with a very low matter density, 
which can hardly be reconciled with the current value derived from the mea- 
surements of the cosmic microwave background anisotropy and galaxy clusters 
(cluster baryon fraction). As a result, this Cardassian expansion proposal does 
not seem to survive the magnitude-redshift test for the present type la supernovae 
data, unless the universe contains primarily baryonic matter. 

Subject headings: cosmology: theory — distance scale — supernovae: type la 
supernovae 



1. Introduction 

A major development in modern cosmology is the discovery of the acceleration of the uni- 
verse through observations of distant type la supernovae (Perlmutter et al. 1998, 1999; Riess 
et al. 1998, 2001; Leibundgut 2001). It is well known that all known types of matter with 
positive pressure generate attractive forces and decelerate the expansion of the universe - 
conventionally, a deceleration factor is always used to describe the status of the universe's ex- 
pansion (Sandage 1988). Given this, the discovery from the high-redshift type la supernovae 
indicates the existence of a new component with fairly negative pressure, which is now gener- 
ally called dark energy, such as a cosmological constant (Weinberg 1989; Carroll et al. 1992; 
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Krauss and Turner 1995; Ostriker and Steinhardt 1995) or an evolving scalar field (referred 
to by some as quintessence) (Ratra and Peebles 1988; Wetterich 1988; Frieman et al. 1995; 
Coble et al. 1997; Caldwell et al. 1998; Gong 2002). While current measurements of the cos- 
mic microwave background anisotropies favor a spatially flat universe with cold dark matter 
(de Bernardis et al. 2000; Lange et al. 2001), both the deuterium abundance measured in 
four high-redshift hydrogen clouds seen in absorption against distant quasars (Buries and 
Tytler 1998a,b) (combined with the baryon fraction in galaxy clusters from X-ray data - see 
White et al. 1993 for the method) and the large-scale structure in the distribution of galaxies 
(Bahcall 2000; Peacock et al. 2001) have made a strong case for a low density universe (for 
a recent summary, see Turner 2002a). It seems that all these observations can be concor- 
dantly explained by the hypothesis that there exists, in addition to cold dark matter, a dark 
energy component with negative pressure in our universe (Turner 2002b). The existence 
of this component has also been independently indicated by other observations such as the 
angular size-redshift relations for compact radio sources (Guivits et al. 1999; Vishwakarma 
2001; Lima and Alcaniz 2002; Chen and Ratra 2003) and FRIIb radio galaxies (Guerra et 
al. 2000; Daly and Guerra 2002; Podariu et al. 2003), the age estimates of old high-redshift 
galaxies (Dunlop et al. 1996; Krauss 1997; Alcaniz and Lima 1999) and gravitational lensing 
(Kochaneck 1996; Chiba and Yoshii 1999; Futamase and Hamana 1999; Jain et al. 2001; 
Dev et al. 2001; Ohyama et al. 2002; Sereno 2002). 

Neither a cosmological constant nor a quintessence, the present candidates for the uni- 
verse acceleration mechanism, however avoid the cosmic coincidence problem - why the 
densities of dark energy and dark matter are comparable today (another related but dis- 
tinct difficulty is the fine-tuning problem, see Carroll et al. 1992 for a discussion of this 
point). Although the tracking field model (Zlatev et al. 1999) provides a possible resolu- 
tion to this problem, a convincing dark energy model with a solid basis in particle physics 
is still far off. Therefore it is desirable to explore alternative possibilities, such as higher 
dimensions (Deffayet et al. 2002; Gu and Hwang 2002) or an altered theory of gravitation 
(Behnke et al. 2002). Very recently, Freese and Lewis (2002) proposed the "Cardassian Ex- 
pansion Scenario" in which the standard Friedman- Robert son- Walker (FRW) equation is 
modified as follows, 

H 2 = Ap + Bp n (1) 

where H = R/R is the Hubble parameter as a function of cosmic time, R is the scale factor of 
the universe and p is the energy density of matter and radiation. In the usual FRW equation 
B = 0. To be consistent with the usual FRW result, one should take A = 8tcG/3. It is 
convenient to use the redshift z eq , at which the two terms of eq.(l) are equal, as the second 
parameter of the Cardassian model. In this parameterization of (n,z eq ), it can be shown 
that (Freese and Lewis 2002), B = F 2 (l + z C9 ) 3(1 " n) Po n [ 1 + (1 + ^ ? ) 3(1_n) ] _1 , where p is the 
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matter density of the universe at the present time and H = 100/?,kms~ 1 Mpc~ 1 is the Hubble 
constant. This particular proposal is very intriguing because the expansion of the universe 
will be accelerated automatically later without any dark energy component - the second 
term, which may arise as a consequence of brane world cosmologies, dominates at a late 
epoch and drives the acceleration of the universe. It is valuable to explore the agreement 
of the Cardassian expansion model with the currently available cosmological observation 
data, as suggested by Freese and Lewis, who proposed this scenario. In a previous paper, 
the authors have used the recent measurements of the angular size of high-redshift compact 
radio sources made by Gurvits et al. (1999) to test the Cardassian model (Zhu and Fujimoto 
2002). It was shown that the allowed intervals for the two model parameters, n and z eq , are 
heavily dependent on the value of the mean projected linear size / (see Table 1 of Zhu and 
Fujimoto 2002). For example, at I = 16h~ 1 pc, the best fit occurs for n = 0.76 and z eq = 1.78, 
which gives a reasonable matter density of Q m ~ 0.32. However, this analysis shows that, if 
one minimizes \ 2 f° r the parameters /, n and z eq simultaneously, the best fit to the current 
angular size data prefers the conventional flat ACDM model to the Cardassian expansion 
proposal. In this work, we analyze this scenario with the distant type la supernovae sample 
compiled by Perlmutter et al. (1999). It is shown that the allowed intervals for n and 
z eq would give rise to a universe with an unreasonably low matter density (Q m < 0.1). 
As a result, the Cardassian expansion scenario does not seem to survive the magnitude- 
redshift test from the present type la supernovae data, unless the universe contains no 
dark matter. After providing a brief description of the data analysis method (section 2), we 
present our numerical results in section 3. Finally we summarize our conclusions and present 
discussion(section 4). 



2. Outline of the data analysis method 

The apparent bolometric magnitude m(z) of a standard candle with absolute bolometric 
magnitude M is related to the luminosity distance di by m = M + 5 logrf^ + 25. Following 
Perlmutter et al. (1997), we write the B-band magnitude-redshift relation as 

m B = M B + 5\ogV L , (2) 

where T>i = H^d^ is the "Hubble-constant-free" luminosity distance and A4b = Mb — 
5 log Hq + 25 is the "Hubble-constant-free" .B-band absolute magnitude at maximum of a 
type la supernova. 

In a Friedmann-Robertson- Walker (FRW) universe, the luminosity distance di can be 
calculated from the redshift-dependent Hubble parameter, H(z) = H E(z), by the integral 



-4- 



cLl{z) = (c/H )(l + z) f* dz'/E(z). For the ansatz of eq.(l) and a flat universe with only 
matter, Freese and Lewis (2002) get 

E\z- n, z eq ) = (1 + (1 + ^) 3(1 - ) )- 1 x (1 + zf + (1 - (1 + (1 + z^f 1 ^)- 1 ) x (1 + zf n (3) 

where n and z eq are the two paramters of the Cardassian model. 

We use the Perlmutter et al. (1999) data to place observational constraint on the 
Cardassian model parameters n and z eq . This data set, plotted in Fig. 1, consists of 42 high- 
redshift Type la supernovae from the Supernova Cosmology Project 1 , and 18 low-redshift 
Type la supernovae from the Calan/Tololo Supernova Survey (Hamuy et al. 1996). Both 
sets include corrections for the lightcurve width-luminosity relation. The error bars, which 
include both the measurement errors and the intrinsic luminosity dispersion, have also been 
shown in Fig. 1. We determine the model parameters n and z eq through a x 2 minimization 
method. The range of n spans the interval [-3, 2] in steps of 0.01, while the range of z eq 
spans the interval [0, 4] also in steps of 0.01. 

2 lKA \ ^[ m B{zi]M B ;n,z eq ) -mfj 2 
X {M B ;n,z eq ) = 2^- 2 ' ( 4 ) 

i a m Bi 

where tob^; Mb', n, z eq ) refers to the theoretical prediction from eq.(2), m|| is the observed 
effective magnitude, and o mm is the total uncertainty (i refers to the ith supernova of the 
sample). The summation is over all of the observational data points. 

Evaluating the ansatz of eq.(l) at the present time, we have (Freese and Lewis 2002) 

^o 2 = ^Po[l + (l + ^) 3{1 ^]. (5) 

Because in the Cardassian model the universe is flat and contains only matter, the matter 
density at present, po, should be equal to the 'critical density' of this scenario. From eq.(5), 
we have 

P0 = Pccardassian = Pc X F(n), F(n) = [1 + (1 + Z eq f {l ~ n) ] ~ L (6) 

where p c = 3Hq/87tG is the critical density of the standard FRW model. Therefore the new 
critical density p c ,cardassian depends on the two parameters, n and z eq , while F(n, z eq ) = Vt m 
gives the matter density in units of the critical density of standard FRW model (Freese 
and Lewis 2002). Instead of specifying Q m (or F), we consider both n and z eq as inde- 
pendent paramters, while Q m (F) is treated as the output of the fitting result. The mag- 
nitude "zero point" A4b can be determined from the 18 low-redshift supernovae that are 
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carefully chosen from a sample of 29 supernovae from the Calan/Tololo survey. After ap- 
propriate correction, they give Mb = —3.32 ± 0.05 (Perlmutter et al. 1997). We will use 
Mb = —3.27, —3.32, —3.37 respectively as typical values of the zero point to fit the data, 
as well as a range of Mb — —3.00 — 3.60 (Meszaros 2002) to check the robustness of our 
results. In order to make the analysis independent of the choice of Mb, we also minimize 
eq.(4) for Mb, n and z eq simutaneously, which we refer as the "best fit". 

3. Numerical results 

Table 1 summarizes our fitting results to the Cardassian expansion model. Following 
Perlmutter et al. (1999), we analyze the 60 supernovae as three different sample groupings. 
Sample A is the entire data set. Sample B excludes four outliers - the two of them with 
lower redshifts, SN1992bo and SN1992bp, are the most significant outliers from the average 
light-curve width, while the other two with higher redshifts, SN1994H and SN19970, are 
the largest residuals from x 2 fitting. Sample C further excludes two very likely reddened 
supernovae, SN1996cg and SN1996cn. (For details of all these outliers, see Perlmutter et 
al. 1999.) As shown in Table 1, the fitting results for samples A and B are very similar 
except for their goodness-of-fits. The larger x 2 P er degree of freedom (d.o.f) for sample A, 
xl = 1.78, indicate that the outlier supernovae included in this sample are probably not part 
of a Gaussian distribution and thus will not be appropriately weighted in a x 2 fit (Perlmutter 
et al. 1999). The x 2 P er d.o.f for sample B, xl — 1-20, is reduced significantly and indicates 
that no large statistical errors remain unaccounted for. The fit for sample C is a more robust 
one, because the two very likely reddened supernovae, SN1996cg and SN1996cn, have been 
further removed (Perlmutter et al. 1999). All three best-fits result in the same value of the 
zero-point magnitude Mb — —3.42, which is higher than Mb = —3.32, the conclusion of 
Perlmutter et al. (1997,1999), but a little bit lower than Mb = —3.45, the value of Efstathiou 
(1999) which is obtained from the best-fit to the combined data for type la supernovae and 
the cosmic microwave background anisotropies. The best fit to sample C, with Mb — —3.42, 
n = —1.33, z eq = 0.43 and the lowest x 2 P er d.o.f. of 1.11, is depicted in Fig. 1 as a solid 
line. For comparison, three other curves with model parameters n and z eq taken from the 
Table 1 of Freese and Lewis (2002) are also shown. The ability of these curves to fit the data 
is surprisingly distinct (and can even be seen by eye), while the former solid curve matches 
the data points very well, none of the later three curves does. 

In Fig. 2, we show the confidence regions (68.3% and 95.4% C.L.) of the fitting results 
in the plane (n,z eq ). The three left panels show the results for sample A, B and C using 
the value of Mb = —3.32 which was initiated by Perlmutter et al. (1997, 1999), while the 
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three right panels show their corresponding best fits (i.e., the case of Mb = —3.42). In 
order to evaluate how reasonable the resulting parameters, n and z eq are, two areas which 
give a currently optimistic matter density, Q m (F) = 0.330 ± 0.035, (see Turner 2002a for 
the argument) and a wider range of Q m (F) = 0.2-0.4, respectively are also shown in every 
panel. As it is distinctly shown that there is no overlap between the resulting parameter 
range and the reasonable area for the matter density, we have a 95.4% confidence level (C.L.) 
in saying that the Cardassian expansion is not compatible with a cold dark matter dominated 
universe with £l m (F) = 0.2-0.4 (the C.L. goes up to 99% if the matter density of the universe 
is, Q m (F) = 0.330±0.035). As a matter of fact, all of our fitting results for n and z eqi point to 
a universe with Q m (F) less than 0.1, which is unreasonable in light of the currently available 
cosmological observations (Bahcall 2000; Peacock et al. 2001; Turner 2002a). 

One could use the observational constraints on the deceleration factor q to cross-check 
the robustness of the fitting results and the difficulty of the Cardassian model with a rea- 
sonable matter density (e.g., Vt m ~ 0.33) in explaining the type la supernovae data. For the 
Cardassian expansion scenario parameterized by n and z eq , we get the deceleration parameter 
as a function of redshift 



where the function E 2 (z;n, z eq ) is given by eq.(3), in which the (1 + ^) 3 -dependent term 
dominates at high redshifts, causing the deceleration of the expansion of the universe. An 
acceleration will only occur when it becomes negative, i.e. q < 0, at a late epoch. We plot 
the redshift dependent deceleration parameter in Fig. 3 for the Cardassian models with the 
parameters of n and z eq taken from the Table 1 of Freese and Lewis (2002) and from our 
best fits. The shaded area in Fig. 3 corresponds to the present observational constraints 
on the deceleration parameter, i.e., the universe switched from deceleration to acceleration 
at a redshift interval 0.6 < z q= o < 1.7 at the la level (Perlmutter et al. 1999; Riess et al. 
1998,2001; Avelino et al. 2001; Avelino and Martins 2002). The problem now is apparent: 
while the best fits of this work predict the turnaround redshift well within the observation 
constraints, all Cardassian models with a reasonable matter density (Q m ~ 0.33) (parameters 
taken from the Table 1 of Freese and Lewis 2002) predict the turnaround redshift less than 
~ 0.6, which is only marginally compatible to the present observations. We are thrown into a 
dilemma: an unrealistically low matter density is needed so that the acceleration starts early 
enough to be realistic (See Avelino and Martins 2002 for the same discussion for another 
brane world cosmology). Therefore our fitting results are robust in explaining the type la 
supernovae data of Perlmutter et al. (1999): all of them predict a universe with very low 
matter density (fi m < 0.1). 
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4. Conclusions and Discussion 

We have analyzed the Cardassian expansion recently proposed by Freese and Lewis 
(2002) using distant type la supernovae data complied by perlmutter et al. (1999). Although 
this particular proposal is an intriguing mechanism for the acceleration of the universe be- 
cause it postulates the universe is flat, matter dominated and accelerating, but contains no 
vacuum contribution, it is strongly disfavored by the present high-redshift type la super- 
novae data and the constraint of Q m ~ 0.3. The main point is that all fitting results of this 
scenario to the supernovae sample lead to a universe with unreasonably low matter density, 
leaving no space for the huge amount of dark matter whose existence has been widely ac- 
cepted among the astronomical community (see, e.g., Primack 2002, Turner 2002b). Even if 
one can say that this Cardassian model can marginally pass the cosmological test from the 
updated angular size data (Zhu and Fujimoto 2002), it can hardly survive the magnitude- 
redshift test for the present type la supernovae data unless the universe contains primarily 
baryonic matter. There seems to be a tendency: a model that excludes the dark energy 
component dispels dark matter also (see Avelino and Martins 2002 for another analysis). 
However, it is worth keeping in mind that a universe with low matter density fl m ~ 0.1 can 
also fit the data of Perlmutter et al. (1999) surprisingly well (Meszaros 2002). 

One of the major uncertainties in the present analysis comes from the errors of the 
magnitude "zero point" Mb- There are several ways to overcome this problem. First of 
all, one can analyzes the data over a large enough range of Mb to include almost all of 
the possibilities, and then calculate the probability distribution for the model parameters by 
integrating over it (Perlmutter et al. 1997, 1999). However this is not needed for our purpose, 
because all our fits with the values of Mb from -3.00 to -3.60 lead to a universe with very low 
matter density. Second, one could pin down the value of Mb through a larger sample low- 
redshift supernovae. Databases of nearby SNela are becoming unprecedentedly abundant 
(see, e.g., Li et al. 2001). In particular, the Nearby SN Factory 2 will accumulate a sample 
of 300 low-redshift supernovae and determine Mb to a precision of ±(0.01-0.02). Finally, 
from eq. (2), the model parameters n and z eq can be determined by measuring differences of 
magnitudes at different redshifts, which are independent of Mb (Frieman et al. 2002). 

Other uncertainties of cosmological parameter extraction from high-redshift type la su- 
pernovae sample caused by progenitor and metallicity evolution, extinction, sample selection 
bias, local perturbations in the expansion rate, gravitational lensing and sample contamina- 
tion have been carefully studied by Riess et al. (1998) and Perlmutter et al. (1999). It was 
found that none of these effects can seriously change the result. Considering our results show 
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that the matter density predicted by the Cardassian scenario is less than 0.1, it is impossible 
for any of these effects to change the case. In short, the present Cardassian expansion model 
is strongly disfavoured by the current distant SNela data and the constraint of Q m ~ 0.3. 
We hope that other convincing mechanisms for the acceleration of the universe will appear 
in the near future. 

We would like to thank A. G. Riess and W. Li for their help on parameter extraction 
with high-redshift type la supernovae data, Y. Tsunesada and S. Sato for their aasistance on 
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help during his stay. Finally, our thanks go to the anonymous referee for valuable comments 
and useful suggestions, which improved this work very much. This work was supported by 
a Grant-in- Aid for Scientific Research on Priority Areas (No. 14047219) from the Ministry of 
Education, Culture, Sports, Science and Technology. 
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Fig. 1. — Hubble diagram for 42 high-redshift SNela from the Supernova Cosmology Project 
and 18 low- redshift SNela from the Calah/Tololo Supernova Survey. The empty squares 
mark the four outliers which are excluded in sample B, while the empty diamonds mark 
the another two futher excluded in sample C. The solid curve corresponds to our best fit to 
sample C, with Mb = —3.42, n = —1.33, z eq = 0.43. The values of (n, z eq ) for the other 
three curves are taken from the Table 1 of Freese and Lewis (2002). 
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Fig. 2. — Confidence regions of fitting results of sample A (a,b), B (c,d) and C (e,f), 
and the area constrained by observed matter density of the universe in the n-z eq plane. 
The faint (dark)-shaded areas show the parameter regions with a confidence level of 68.3% 
(95.4%). The positively slanted hatchings and the cross-hatched regions correspond to the 
parameter areas that give the matter density of the universe Q m (F)— 0.2- 0.4 and Q m 
(F)— 0.330 ± 0.035 respectively. The left three panels are for the fitting results with the 
assumption of Mb = —3.32, while the right three panels show the case of Mb = —3.42. 
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redshift z 

Fig. 3. — Diagram of the deceleration parameter q versus redshift for the Cardassian expan- 
sion models. The model parameters are taken from Table 1 of Freese and Lewis (2002) and 
from our best fits respectively. The shaded area shows the observational constraint on q from 
the literatures. As it shows, our best-fits are much more compatible with the observational 
constraints on the turnaround redshifts than the Cardassian models with Q m (F) = 0.33 are. 
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Table 1. Fitting results for the Cardassian model from distant type la supernovae data 

compiled by Perlmutter et al. (1999). 
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0.49 


0.068 


102.1 


A 


60 


Best fit: 


-3.42 


-1.43 


0.39 


0.083 


101.3 


B 


56 




-3.27 


-0.80 


0.83 


0.037 


69.7 


B 


56 




-3.32 


-0.99 


0.61 


0.055 


66.4 


B 


56 




-3.37 


-1.21 


0.47 


0.072 


64.5 


B 


56 


Best fit: 


-3.42 


-1.45 


0.38 


0.086 


63.9 


C 


54 




-3.27 


-0.71 


1.08 


0.023 


65.4 


C 


54 




-3.32 


-0.90 


0.72 


0.043 


62.4 


c 


54 




-3.37 


-1.12 


0.53 


0.063 


60.7 


c 


54 


Best fit: 


-3.42 


-1.33 


0.43 


0.076 


60.2 



Note. Sample A: all supernovae; Sample B: excludes 
four outliers, SN1992bo, SN1992bp, SN1994H and SN19970; 
Sample C: further excludes SN1996cg and SN1996cn. 



